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Abstract. Perhaps the most controversial idea in modern cosmology is that our observable 
universe is contained within one bubble among many, all inhabiting the eternally inflat¬ 
ing multiverse. One of the few way to test this idea is to look for evidence of the relic 
inhomogeneities left by the collisions between other bubbles and our own. Such relic inho¬ 
mogeneities will induce a coherent bulk flow over Gpc scales. Therefore, bubble collisions 
leave unique imprints in the cosmic microwave background (CMB) through the kinetic Sun¬ 
yaev Zel’dovich (kSZ) effect, temperature anisotropies induced by the scattering of photons 
from coherently moving free electrons in the diffuse intergalactic medium. The kSZ signature 
produced by bubble collisions has a unique directional dependence and is tightly correlated 
with the galaxy distribution; it can therefore be distinguished from other contributions to 
the CMB anisotropies. An important advantage of the kSZ signature is that it peaks on 
arcminute angular scales, where the limiting factors in making a detection are instrumental 
noise and foreground subtraction. This is in contrast to the collision signature in the primary 
CMB, which peaks on angular scales much larger than one degree, and whose detection is 
therefore limited by cosmic variance. In this paper, we examine the prospects for probing 
the inhomogeneities left by bubble collisions using the kSZ effect. We provide a forecast 
for detection using cross-correlations between CMB and galaxy surveys, finding that the de¬ 
tectability using the kSZ effect can be competitive with constraints from CMB temperature 
and polarization data. 
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1 Introduction 

The idea that we might inhabit an eternally inflating multiverse is surely one of the most 
radical concepts to come out of modern physics. Although radical, accelerated expansion 
of the universe and spontaneous symmetry breaking, both observed phenomena in nature, 
are the necessary and sufficient ingredients. Once accelerated expansion driven by vacuum 
energy begins, an epoch known as inflation in the early universe, it is generically difficult 
to end globally. Instead, inflation ends only locally inside of bubbles which are formed only 
rarely via quantum mechanical effects. If bubbles are formed slower than the space between 
them is expanding, inflation cannot end globally, giving rise to the phenomenon of eternal 
inflation. 

In this picture, our universe is housed within one bubble universe among many, perhaps 
each with different properties and cosmological histories (as motivated by the string theory 
landscape [1]). The Coleman de Luccia instanton [2-4] that mediates bubble nucleation 
dictates the symmetry of a one-bubble spacetime, giving rise to an open FRW universe in 
the bubble interior. If we live inside a bubble, a second epoch of slow-roll inflation is necessary 
to dilute curvature to observable levels and provide the seeds for structure formation [5-7]. 

Each bubble undergoes collisions with others [8] , producing potentially observable wreck¬ 
age [9] . Assessing the character of this wreckage and the probability of observing it has been 
the subject of a substantial body of work [10-43]. Given an underlying model consisting of 
a scalar field theory coupled to gravity, a full set of predictions can be made for the primor¬ 
dial metric fluctuations caused by a bubble collision [25, 29]. The symmetry of the collision 
spacetime forbids the production of gravitational waves [11, 12], and therefore the effects 
of a collision in any single-field model can be encoded in the primordial scalar comoving 
curvature perturbation. Bubble collisions can be mapped onto specific inhomogeneous initial 
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Figure 1. The scalar field potential underlying the bubble collision model. 


conditions for slow-roll inflation inside of the bubble. Because the wreckage produced by 
bubble collisions is pre-inflationary, having the minimum number of inflationary e—folds is 
necessary for observable effects to be accessible today 

Each collision is an event, and therefore can only affect regions in its causal future. The 
bubble interior is split into regions that are affected by the collision and regions that are 
not. Observers in the vicinity of this causal boundary have access to the most distinctive 
observational signatures; in the following we will focus on this class of observers. The first 
constraints on cosmic bubble collisions in the eternal inflation scenario were found by Feeney 
et al [46-49] using observations of the Cosmic Microwave Background (CMB) radiation made 
by the WMAP satellite [50]. This was followed up by the analysis of Osbourne et al [51, 52]. 
Neither of these analyses support the bubble collision hypothesis, placing a constraint on 
the possible models of eternal inflation. However, these analyses were performed with an 
incomplete understanding of the collision signature and the connection between the under¬ 
lying theory and observable predictions. An analysis incorporating the improved predictions 
of [25, 29] is underway. 

Because the dominant contribution to the CMB temperature anisotropies arises from 
the time of last scattering, it is not an ideal probe of inhomogeneity. CMB polarization probes 
the local temperature quadrupole observed by electrons, providing important additional in¬ 
formation on inhomogeneities (see e.g. [53, 54]) and acting as an important discriminator 
between anisotropy and inhomogeneity. However, since the signature of a bubble collision 
is imprinted on the largest scales, cosmic variance is a significant impediment to placing 
stringent constraints on the bubble collision model. Can we ever hope to do better? 

This paper explores the ability of another probe of large-scale inhomogeneities, the 
kinetic Sunyaev Zel’dovich (kSZ) effect [55], to test the bubble collision hypothesis. The kSZ 
effect is caused by the inverse Compton scattering of CMB photons by coherently moving 
free electrons. Its power to probe missing baryons, reionization and peculiar velocity at 
< 100 Mpc scale has been long recognized (e.g. [56-60]), and has been explored in CMB 
experiments such as ACT and SPT (e.g. [61, 62]). It can also be utilized to constrain dark 
matter-dark energy interaction [63] and modified gravity [64, 65]. Recently, a new regime 
of the kSZ effect generated by larger scale (> Gpc) inhomogeneities has been realized and 
applied to probe dark flows [66-72], test the Copernican principle [73, 74], and constrain the 

^Various arguments may be advanced for the plausibility of this situation, see e.g. [44, 45] 
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Figure 2. A schematic diagram of the kSZ effect in the collision spacetime. The causal boundary of 
the collision Xc is denoted by the dotted line (plane). The solid circle (sphere) is the last scattering 
surface (LSS) seen by us, and the dashed circles (spheres) are the LSS of several representative free 
electrons, labelled by “el”, “e2” and “e3” respectively. These electrons see corresponding CMB 
dipoles on their LSS and have their own peculiar motion with respect to our CMB. In contrast to 
a single adiabatic superhorizon mode, the bubble collision generated perturbation induces ?;eff 0 
and therefore induces a non-vanishing kSZ effect. Unlike the primary CMB fluctuation generated 
by bubble collision, which cover only a fraction of the sky, the kSZ effect generated by the bubble 
collision is spread over the whole 47r sky. For example, all electrons at a radius between “e3” and us 
(“O”) at all viewing angles can see the bubble collision through their LSS and hence generate a kSZ 
effect visible to us. 


present day vacuum decay rate [75]. This paper presents its newest application. 

The kSZ effect arises from scattering of CMB photons from free-electrons undergoing 
bulk motion in the late universe. The amplitude of the effect is proportional to the locally 
observed CMB dipole, and is therefore a sensitive measure of bulk velocities.^ Horizon 
scale inhomogeneities generate bulk motions coherent over Gpc scales with characteristic 
directional dependence, and therefore, a distinctive kSZ signal. Note that this is a special 
class of kSZ signatures, distinctive from the conventional kSZ effect used to probe peculiar 
velocities on small scales. 

The kSZ effect is a unique and powerful probe of horizon scale inhomogeneities, for a 
number of reasons. 

1. Since the distance to the epoch of reionization (at z ~ 6—10) is very close to the distance 
to our horizon, free electrons probe much of the observable universe. Therefore the 
kSZ effect provides a unique opportunity to observe the universe from virtually every 

^Although we do not consider it in this paper, polarization of the kSZ photons map the locally observed 
CMB qnadrnpole [76]. 
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position inside of our horizon, and is sensitive to horizon scale inhomogeneities such as 
those generated by a bubble collision. 

2. The kSZ signature is essentially a modulated re-mapping of the electron distribution to 
the CMB. Crucially, due to small scale inhomogeneities in the electron distribution, the 
power of the kSZ effect generated by horizon scale inhomogeneities peaks at arcminute 
angular scales.^ Hence, the kSZ effect promises to give constraints beyond the cosmic 
variance limited results obtained from the primary CMB only. 

3. In contrast to the conventional kSZ effect, whose two-point cross correlation with the 
galaxy density field vanishes, the new kSZ effect is tightly correlated with the galaxy 
density field [67, 74]. Furthermore, the overall amplitude of the cross correlation has a 
unique directional modulation, depending on specific models of horizon scale inhomo¬ 
geneities. Hence by specially designed cross correlation with galaxies, known as kSZ 
tomography [77], we can not only isolate it from other CMB components (primary 
CMB, other kSZ components, cosmic infrared background, etc), but also determine its 
redshift dependence. 

The remainder of the paper is structured as follows. In Sec. 2 we discuss the predicted 
inhomogeneities produced by bubble collisions. Sec. 3 computes the kSZ signature expected 
for a general bubble collision. Sec. 4 estimates the capability of probing bubble collision using 
kSZ tomography, and makes a comparisong with constraints obtained from the primary CMB. 
Finally, we discuss and summarize in Sec. 5. 

2 Cosmology in the aftermath of a bubble collision 

A model of eternal inflation is defined by a scalar field lagrangian, which by assumption has 
a metastable false vacuum and at least one true vacuum. Transitions from the false to the 
true vacuum are governed by the CDL instanton, which gives rise to bubbles of true vacuum 
nucleating in the false vacuum. The cosmological history inside the bubbles includes an epoch 
of cosmic inflation, and the 0(4) symmetry of the CDL instanton guarantees that a single 
bubble contains an infinite open Friedman Robertson Walker (FRW) universe. Focusing 
on a single bubble, which we hypothesize might contain our observable universe, additional 
bubbles forming from the false vacuum will collide with ours forming a fractal cluster of 
intersecting bubbles [8, 78]. In the limit where collisions are rare, they will lie outside of 
each other’s horizon; we can therefore treat collisions on a pair-wise basis. The spacetime 
describing the collision between two bubbles has an S0(2,l) symmetry [12], inherited from 
the symmetry of each of the colliding bubbles. Determining the precise outcome of a collision, 
and therefore the observable effects, has been examined using analytic and numerical meth¬ 
ods. Recent efforts by Wainwright et al [25, 29] have established a connection between the 
primordial curvature perturbation and specific lagrangians, forming a first-principles deriva¬ 
tion of the effects of bubble collisions in single field models. The conclusions of this work are 
as follows: 

^Roughly speaking, fluctuation in such kSZ effect is proportional to the matter density fluctuation projected 
along some radial depth. Therefore the peak scale depends on the linear transfer function, ag, relevant 
redshifts, and the projection length. For typical LCDM density fluctuation projected along Gpc radial depth, 
the typical peak scale is of the order one arcminute. This is numerically demonstrated in cases of dark flow 
induced kSZ effect [67]) and the kSZ effect discussed in this paper (Fig. 9). 
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• The boundary of the causal future of the collision splits the inner-bubble cosmology 
into two qualitatively different regions. This work considered observers who are in the 
vicinity of this boundary, reconstructing the perturbed FRW universe in this region in 
comoving gauge. 


• On the scale of our observable universe, where spatial curvature can be neglected, the 
collision spacetime has approximate planar symmetry. 


• Neglecting spatial curvature (at most a 1% effect), the perturbed FRW universe in 
Newtonian gauge is 

ds"^ = -(1 + 2^)dt‘^ + - 2T)dx2 (2.1) 


with 


^i(r) ~ 


A{x — Xc) + B{x — XcY, if X > Xc 
0, if X < Xc . 


( 2 . 2 ) 


Here, we have fixed the collision direction to be along the x axis, and Xc is the location 
of the causal boundary. The prior probability distribution over Xc is flat over the size 
of the present day horizon, x = x/vh and rn = is the present day Hubble radius 

(where c is the speed of light). 


• A comparison between numerical and analytic methods in [29] showed that the coeffi¬ 
cients A and B can be related to fundamental parameters of the scalar field lagrangian 
shown in Fig. 1 by 


A = -1 


robs Mpi 


(1 - cos AXsep) 


(2.3) 


B = 


^_Qobs /'^coll Hj 


15 


coll 


^’obs Hf 


obs 


(1 - cos AXsep)" 


where 64>o is the distance between the instanton endpoints connecting the false vacuum 
to the collision bubble interior, is the energy density in curvature in our own 

universe (constrained to be < 10“^), Xobs is the tensor to scalar ratio for inflation 
inside the observation bubble (constrained to be robs ^ -1), Xcoii is the tensor to scalar 
ratio for inflation inside the collision bubble, is the Hubble scale during inflation 
inside the observation bubble, is the Hubble scale during inflation inside the 

collision bubble, and Axgep is the initial proper distance between the bubbles in the 
collision frame measured in terms of the false vacuum Hubble scale (0 < Axsep < ^r)- 
By definition, 64>q is positive for collisions between identical bubbles and negative for 
collisions between different bubbles. The prior probability distribution over initial 
separations Axgep is proportional to sin^ Axgep- We therefore expect (1 — cos Axgep) ~ 
0(1). The analytic expression Eq. 2.3 assumes that slow-roll inflation is not disrupted 
by the collision, in which case is roughly constant everywhere in the collision 

spacetime during inflation. Complete disruption of inflation will occur when is 

comparable to the total field excursion during inflation inside the bubble. Using the 
Lyth Bound [? ] (yT^ = y/0.0lA(j)/Mpi), we therefore conclude that Eq. 2.3 will be 
valid as long as < Mpiy^robs/O.Ol. 


■^Comparing with Ref. [29], a factor of 3/5 has been applied to convert between the Newtonian and comoving 
gauges, and we have converted our dimensionless measure of distance from the distance to the surface of last 
scattering to the size of the present day Hubble radius, introducing another factor of 3. 
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• We restrict ourselves to so-called native born observers [16], that is those who are 
comoving with respect to the unperturbed portion of the bubble. By definition, this 
restricts Xc > 0. There are also observers comoving with respect to the perturbed 
portion of the bubble. In the vicinity of Xc, these observers would have causal access to 
the collision boundary. Here, and also deep within the collision region, the dominant 
perturbation would be a quadratic curvature perturbation centred on the observer. 

Cosmological datasets can be used to constrain the parameters in the template Eq. 2.2, 
which in turn constrains the underlying scalar field lagrangian through Eq. 2.3. We now turn 
to assessing the ability of the kSZ signature produced by bubble collisions to be used as such 
a probe. 

3 The kinetic Sunyaev Zel’dovich effect generated by bubble collision 

The curvature perturbation induced by bubble collisions gives rise to relative motion between 
matter in different regions of the Universe. Inverse Compton scattering off free electrons in 
relative motion induces temperature anisotropies in the observed CMB sky via the kSZ effect. 
The relevant induced temperature anisotropies are given by [55] 

(he) = -or f dXe Oe ne{x{Ze)n, Ze) . (3.1) 

J-CMB Jo C 

The minus sign is chosen such that when an electron moves away from us (u > 0), the kSZ 
effect is negative. The vector hg points from us to the CMB sky, which is also the vector from 
us to the free electrons scattering CMB photons. The subscript “e” denotes properties of 
free electrons. Og = l/(l-|-^;g) is the scale factor. Xg is the comoving radial coordinate of free 
electrons, and dxe = cdzejH. For a flat universe, Xg coincides with the comoving angular 
diameter distance, ng is the electron number density. Ugff is the effective velocity generating 
the kSZ effect, which we derive in detail below. Besides the peculiar velocity of electrons 
at redshift z and direction hg, it also takes into account the contribution from the intrinsic 
CMB dipole on the last scattering surface (LSS) of free electrons at different locations. This 
new contribution only exists in non-Copernican universes. 

A schematic plot of the kSZ effect produced in the bubble collision spacetime is shown 
in Fig. 2. Since the distance to the reionization epoch at z ~ 10 is very close to the distance 
to our horizon, free electrons are essentially everywhere inside of our horizon. Each electron 
sees its own LSS, and combining all these LSS enables us to fairly sample all regions inside 
of our horizon. For example, the primary CMB along the line of sight from “O” to “e3” 
is not affected by the depicted bubble collision. Nevertheless, the kSZ effect generated by 
electrons located between “e3” and “O” allows us to infer the existence of a bubble collision 
along this line of sight. At fixed distance from us, the kSZ effect is confined within a disk. 
However, combining electrons from all distances spreads the effects of the bubble collision 
over the whole sky. 

3.1 The effective velocity responsible for the kSZ effect 

To calculate Ugff, we need to determine the evolution of the Newtonian potential and velocity. 
The perturbation induced by bubble collisions (Eq. 2.2) is a mixture of both superhorizon 
and subhorizon modes, with a dominant contribution from superhorizon modes. Subhorizon 
modes are further damped compared to superhorizon modes during cosmological evolution. 
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with the transfer function satisfying T[k ;§> 1/H)/T[k <C ^/H) ^ 1- Therefore we expect 
that the overall evolution of the potential and peculiar velocity can be approximated by that 
of the superhorizon modes, which are scale independent, and have analytical expressions (e.g. 
[79]). We outline this computation in Appendix A. To test this approximation, we numer¬ 
ically compute the potential and peculiar velocity by solving the Boltzmann equations for 
radiation and cold dark matter (CDM) in a flat universe. We confirm that the superhorizon 
approximation works very well (Fig. 17, Appendix A). Including the sub-horizon evolution 
acts to slightly smooth the causal boundary at Xc, a consequence of the damping of sub¬ 
horizon (small scale) modes. We neglect this effect, as well as the baryon-photon coupling 
before last scattering Hereafter we will adopt the super horizon approximation and use the 
analytic expressions in Appendix A to calculate the kSZ effect induced by bubble collision. 

The potential T at late times follows the same form as Eq. 2.2 but with a scale-invariant 
growth factor Dx^{a) (given by Eq. A.6): 

T(r,a) = T),i,(a)Tj(r) . (3.2) 


Eor the velocity field, only the x-component is non-zero. It is well described by 

, . / N -1 \ A + 2B{x — Xc), if X > Xc 

V:,{v,a) = -D^{a)rjlx{ \ “ 

" I 0, if X < Xc . 


(3.3) 


The linear velocity growth factor Dy is given by Eq. A.8. Throughout the paper we will 
discuss two limiting cases with {A,B) = (/ 0,0) (case A) and {A, B) = (0,/ 0) (case B). 

The effective velocity in Eq. 3.1 is related to the CMB dipole observed by a freely falling 
electron via 

Vesihe) = ^ d^n 0(rdec) (n • fie) ■ (3.4) 

Here cos 9 = h-hg and the unit vector h points from the electrons to their LSS. r^ec = Xehe + 
^Xdech. Here Xe = xi^e) is the distance from us to free electrons, and Aydec = x(^dec) — Xe 
is the distance from the electrons to their LSS. The geometry is shown in Fig. 2. The effective 
velocity has three contributions. 


Vesifle) = ^^eff,SW + Doppler + t’eff,ISW • (3-5) 

These contributions come from the Sachs-Wolfe (SW) effect generated by the gravitational 
potential on the LSS, the Doppler effect due to peculiar motion of photons on the LSS and 
peculiar motion of electrons at redshift Ze, and the integrated Sachs-Wolfe (ISW) effect [80], 
respectively. For a freely falling electron at redshift Zg and position rg = Xef^e, the CMB 
temperature it sees along the direction n is 


©(n) = 0sw(h) + ©Doppler (h) -|- 0ISw(h) . 


(3.6) 


We follow [79] to evaluate the above terms. The SW term is 


0sw(h) 


2 


5 9/10 

3 D<j>(udec) 


'h(''dec, ®dec) 


2 


5 9/10 

3 T)$(Q,dec) 


i4\]/(ffldec)'hi(j^dec) • 


(3.7) 


®This gives rise to the “cosmic wakes” described in Ref. [23], and may add structnre to baryon peculiar 
velocities on BAO scales in the vicinity of the causal boundary beyond what is considered here. We do not 
expect such signatures to significantly alter our results. 
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Figure 3. The effective peculiar velocity that generates the kSZ effect, for Zc = 0. Electrons sit 
at Ze = 0,1,2,4 respectively. The solid lines have B = 1 and A = 0. The dash lines have A = 1 
and B = 0. Zc is the closest redshift to the causal boundary of the collision. Due to the azimuthal 
symmetry of the initial perturbation, reS = 0 when 9^ = tt/2. For A = 1 and i? = 0, we have the 
analytical expression at Zg = 0: VeS = VeCOs{9e) when 6^ > 7r/2 and —VeCOs(9e) when 9e < 7r/2. 
The value of Ve is given in Fig. 17. For Ze = 4, VeS approaches zero when 0e 0 or ^ tt, but 
for different reasons. The reason it approaches zero when ^ tt is that the electrons and their last 
scattering spheres are located in unperturbed regions. The reason it approaches zero when ^ 0 is 
due to the cancellation of SW, ISW and Doppler components [79] . 


Here we have neglected the impact of 'h(re,ae) since it does not cause an observable effect. 
The Doppler effect has two contributions with opposite sign, 

0Doppler(fi) = n ■ [v(re,ae) - v(rdec, Odec)] = n ■ (Odec)VTj(fdec) “ D^(Oe)VTi(re)] (3.8) 

The ISW effect is given by 

— (r(a),a)da = 2 / -^x[>,(r(a))da . (3.9) 

dec ^dec 

Here, r(o) = Fg + A^h. Ax(a) = da/a^H. So Axdec = Ax(adec)- 

Before proceeding, there is an important issue to clarify. For adiabatic perturbations in a 
flat universe, and in the limit of constant gradient in the gravitational potential corresponding 
to a very superhorizon mode, the three contributions (SW, Doppler and ISW) cancel exactly 
(0(h) = 0) [79, 81-83]. Therefore Ueff = 0 (vanishing CMB dipole). This was first found 
by [81-83] for a Dm = 1 Einstein-de Sitter universe. Erickcek et al. [79] proved that it 
also holds in an ACDM universe with radiation and an CDM universe with smooth dark 
energy and radiation. They further argue that this should be a generic feature for adiabatic 
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Figure 4. Same as Fig. 3, but for Zc = 1- Notice that when 9^ ^ tt and > 2, the electrons 
and their last scattering spheres are located in unperturbed regions and hence Ueff —>■ 0. The sharp 
decrease at 9e ~ 50° for Ze = 2 (and at de — 60° for Ze = 4) and B = 0 is caused by the discontinuity 
of velocity at the edge of bubble collision. This sharp feature shows up when Zg > Zc- Ves{n) shows a 
rich dependence on the observation angle, a precious property for clean extraction from an otherwise 
noisy CMB sky. 


superhorizon modes so that the cancellation holds in all flat cosmologies. For these cases, 
the kSZ effect vanishes. Fortunately for the bubble collision generated inhomogeneities, the 
exact cancellation no longer works, for two reasons. First, the curvature perturbation has a 
quadratic term, and hence a spatially varying gradient. So in general Ueg / 0. Second, even 
if the quadratic term vanishes {B = 0), due to the discontinuity in gradient on the edge of 
a bubble collision, Ueff 7 ^ 0 for electrons whose LSSs are not completely Inside of the region 
perturbed by bubble collision (e.g. at position “e3” in Fig. 2). Although electrons whose LSS 
are completely inside of the region affected by bubble collision (e.g. “el” and “e2” in Fig. 2) 
can have Ueg = 0 for B = 0, overall Veg 7 ^ 0 and we expect a non-vanishing kSZ effect. 

The azimuthal symmetry of the collision spacetime allows us to choose the parameteri¬ 
zation 


he = (cos 06) sin^e, 0) , (3.10) 

h = (cos 0 ,sin 0 cos(/?, sin 0 sin(/ 9 ) . 

The integral over y? in Eq. 3.4 is always [0, 27r]. To evaluate the integral over 6 in Eq. 3.4, we 
must consider three cases: 

1. Xe cos 0e + Axdec < Xc- The LSS of the electron is outside of the region affected by the 
bubble collision and therefore Ueff = 0 . 
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Figure 5. Same as Fig. 3, but for Zc = 2. Notice that, for Ze = 0, v^s is anti-symmetric under the 
operation o tt — 0e, a direct result of azimuthal symmetry. The only exception is when = 0 
and 13 = 0 where the anti-symmetry is reversed to symmetry, caused by the sharp transition in the 
velocity field at Ze = Zc = 0 (Fig. 3). 


2 . XeCosBe — Axdec > ^c- The LSS of the electron is completely inside of the region 
affected by the bubble collision. So the integral is over B G [0, tt], 

3. XeCosOe — Axdec < Xc & Xe COS 0e + Axdec > Xc- Only part of the LSS is inside of the 
region affected by the bubble collision. Therefore the integral is over [0,0c]- &c is the 
solution of 9 to Xe cos 9^ + Axdec cos 9 = Xc-, namely cos 9c = {xc — Xe cos 9e) /Axdec- 

For cases 2 and 3, we obtain 


^^eff,SW = 


2-5 


9/10 


3 (ffldec) 


Ti\I/(ndec) 2 COS 0e 


... ^ ,COS^0 ^ COS^0 

A (Xe COS 9e - Xc) -r-h Axdec- 5 - 


7-)// /I .oCOS^^ , „ .. 2cos^0 AxLrCOS^0 

+ -B (XeCOS0e - Xc) — -h (Xe COS 0e - Xc)Axdec- 5 -^-- 


0 

(3.11) 

0 c 


^eff,Doppler 


= v(re, Oe) • ’le + 2 COS 9eDy{adec) 

A COS^ 9 ^ 2B{Xe cos 9c — Xc) COS^ 9 ^ BAxdec COS'^ 9 


(3.12) 
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Figure 6. Simiar to Fig. 3, but for Zc = 4 and Zg = 0, 2,4,6 instead. The two cases of initial 
fluctuations produce similar Ueff at ^ Zc^ but this is a coincidence. At higher Zg, the two cases 
differ significantly. In particular, the sharp transition at Ze = 6 and B = 0 shown in Fig. 5 shows up 
here when Ze = 6 . These electrons sit at the edge of bubble collision, satisfying Xe cos^e = Xc = xi^c)- 


To calculate the ISW effect, we just need to replace Tj in Eq. 3.9 by 
3 cos 6e 




... ^ , cos^ 6 cos^ 9 

A (Xe cos 9e - Xc) — -h Ax „ 


^ .2 008^9 2cos^0 Ax^cos‘^0 

+ (Xe cos 9e - Xcy— -h (Xe COS 9^ - Xc)Ax -r-^- - - 


0 

0c 


Notice that now cos 9c = {xc — Xe cos 9c)/Ax{a). In the limit of no intrinsic CMB dipole, we 
recover the conventional kSZ effect with Ueg = v(re,ae) • he- 

Fig. 3, 4, 5, 6 , 7 & 8 show the values of Ueff and its dependence on 9e, Zg, Zc, A and 
B. Zc is the minimum redshift to the edge of bubble collision, with Xc = xi^c)- For quick 
reference, we plot the Xc-Zc relation in Fig. 16. As expected, Vce/c is in general of the same 
order as A and B. Nevertheless, numerical results show that the maximum value of Vcs/c 
is often smaller than A or i? by a factor of a few to 10. This is likely caused by incomplete 
cancellation between the SW, ISW and the Doppler effect [79, 81-83]. 

We show Ueff as a function of 9c for different values of Zg = 0, 1, 2 &: 4, and Zc = 0 
(Fig. 3), Zc = 1 (Fig. 4), Zc = 2 (Fig. 5) and Zc = 4 (Fig. 6 ). Ueg shows unique 
directional dependence on 9c, which in general can not be mimicked by conventional sources 
and contaminants. Furthermore, such directional dependence relies on Zg, Zc, A and B in 
complicated and different ways. For the purpose of extracting the bubble collision generated 
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Figure 7. The VeS-Ze relation for 4 . = 1, = 0 and electrons with 0^=0 (solid lines), 7r/4 (dot 

lines), 7r/2 (dash lines), 37r/4 (long dash lines), and tt (dot-short dash lines), respectively. Notice 
that due to the azimuthal symmetry, Veff(0e = 7r/2) = 0. For Zc = 0, the lines with 0^=0 overlap 
with those for 0e = tt. This can be understood by the result that tigff = 0 if = A(x — Xc) without 
the causal boundary a.t x = Xc [79]. Also for this reason, the lines 0^ = 7r/4 overlap with those for 
0e = 37r/4. An interesting behavior is that for 0f, = 0, Ve(f(ze) is nearly a constant when Ze < Zc- 


kSZ effect from otherwise overwhelming contaminants and constraining collision parameters 
(4, B and Xc = xi^c)), this is good news. 

Although in general these features can only be obtained numerically, some of them can 
be understood analytically. 

1. For example, due to the azimuthal symmetry, Ueff(0e = '^/2) = 0. This holds for all 
values of Ze, Zc, A, B and Xc- Also for this symmetry, Uefr('/r — 6e) = —Ves{9e) when 
Ze = 0 and Zc > 0. When Zc = 0, this symmetry still holds for the case B. But for case 
A, due to the discontinuity in the velocity field at Zg = 0, we have Ueff (tt — 0e) = Ves{9e) 
instead. 

2. For case A with Zc = 0 and Zg = 0, we have Ves{9e) = Ve{z = 0) cos Be when Be > vr/2, 
and Ves{Be) = —Ve{z = 0) cos Be when Be < vr/2. The value of Ve{z = 0) is given in Fig. 
17. 

3. For some combinations of Zg, Zg and 0g, the LSS of electrons is located outside of the 
region affected by the bubble collision. Therefore we have Ugfr = 0 there. 

4. For some other combinations of Zg, 2 g and 0g, the LSS of electrons are located completely 
inside of the region affected by the bubble collision. Therefore we have Ves = 0 for case 
A, due to the cancellation of the SW, Doppler and ISW effect [79, 81-83]. 
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Figure 8. Similar to Fig. 7, but for the VeS-Ze relation with B = 1 and A = 0. 


5. For case A, there is often a discontinuity in — 0^. This happens when electrons with 
varying Oe cross the causal boundary of collision, where the velocity is discontinuous. 

It is also useful to examine the dependence of Ueg on 2e for fixed Og. This, together with 
ahe oc (1 + z)^ and the evolution of electron number overdensity, determines the contribution 
of electrons at a given redshift to the kSZ effect. Therefore we plot the v^s-Ze relation in 
Fig. 7 & 8, for various Og = 0,7r/4,7r/2, Svr/d &: tt. As expected, Ugfr along a fixed line of 
sight is coherent over Gpc scales. It thus avoids cancellation of velocity of different signs in 
the conventional kinetic SZ effect [56]. Also for this reason, the kSZ-galaxy number density 
cross correlation does not vanish, in contrast to the vanishing cross correlation between the 
conventional kSZ and galaxies (e.g. Fig. 1, [77]). 

3.2 The kSZ auto power spectrum 

The kSZ auto power spectrum induced by a bubble collision is extremely small. For an order 
of magnitude estimate, we can compute 

AT/T ~ T{Vefi/c)6e (3.13) 

with r ~ 0.01 and density fluctuation Jg ~ 0.1 at arcminute scales where the signal will peak. 
From the previous section, the maximum value of Ueff/c is roughly set by the magnitude of 
A and/or B. As we discuss in more detail in Sec. 4.3, existing constraints from the collision 
signature in the primary CMB are roughly A < 10“^ {B has yet to be constrained in the 
primary CMB). The level of the signature in kSZ power spectrum is therefore AT^ < 1/iK^ 
(for the unrealistic values A,B^ 10“^) to AT^ < .Ol/xK^ (at the upper bounds allowed by 
the primary CMB). This is far smaller than the most recent measurement of the kSZ power 
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Figure 9. The kSZ power spectrum for B = 10“^ along the line of sight towards the centre of 
the bubble collision. Solid, dot, shot dash and long dash lines have Zc = 0,1, 2,4 respectively. The 
generated kSZ effect is sub-dominant to the conventional kSZ effect calculated using the model of [59]. 
It is also well below the measured band-power of the conventional kSZ effect at £ ^ 3000 (2000 < 
I < 10000) by SPT [62]. Given the existence of overwhelming contamination from a combination of 
the primary CMB, conventional kSZ and cosmic infrared background, it is infeasible to detect bubble 
collision through a measurement of the kSZ auto power spectrum. 


spectrum from the South Pole Telescope [84] (SPT) which finds the bandpower AT^ = 
2.9 ± 1.3/rK^ Sit I = 3000, notwithstanding the primary CMB, thermal SZ, cosmic infrared 
background, and instrumental noise. 

To confirm this rough estimate, we can perform a more accurate computation following 
the derivation of the dark flow induced kSZ effect outlined in Ref. [67]. Using the Limber 
approximation, the kSZ power spectrum in the he direction is given by. 






z) 

'VesincZeY 

VXe’ V 

lOkm/s 




+ ZefXedXe ■ (3.14) 


Ag(/c,z) is the power spectrum (variance) of the electron number overdensity at redshift 
z. In principle it has contributions both from the conventional primordial fluctuation and 
from the bubble collision. Observationally we know that the first dominates otherwise the 
concordance ACDM cosmology would not work so well. Therefore we will approximate Ag 
as that predicted by the standard ACDM cosmology.® 

Fig. 9 shows the numerical results for the kSZ auto power spectrum for case B with 
B = 10“^ toward the center of the bubble collision. It peaks at ^ ~ 10^, with peak amplitude 
< l^K^. It is subdominant to the conventional kSZ effect by a factor of a few at all angular 

®We have numerically verified that this is indeed the case for interesting values ot A, B < 10“®. 
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Figure 10. The kSZ power spectrum for A = 10 Solid, dot, shot dash and long dash lines have 
9f, = 0°,45°, 135°, 180° respectively. Notice that when 0e = 90°, the kSZ effect vanishes. 


scales. It is also subdominant to the lensed primary CMB at £ < 4000. Fig. 10 & 11 
show that this sub-dominance holds for other lines of sight and for both case A and B, 
with A,B = 10“^ respectively. The auto power spectrum for other values of A, B scales as 
(A, B/10~^)‘^. Since only A,B< W~‘^ is allowed by constraints from primary CMB, the kSZ 
effect induced by bubble collisions cannot be detected in the temperature power spectrum. 
Nevertheless, we now show that kSZ tomography, namely the kSZ-large scale structure cross 
correlation, is a promising probe for bubble collisions. 

3.3 The kSZ-large scale structure cross correlation 

Although the kSZ signature of bubble collisions in the temperature power spectrum is hope¬ 
lessly buried, kSZ tomography can be used to boost the signal by specifically designed 
weighted cross correlations of CMB temperature with galaxies [67, 77, 85, 86]. A schematic 
is shown in Fig. 12; in a set of redshift bins the kSZ effect caused by bubble collisions is a 
modulated map of the galaxy distribution onto CMB temperature. This technique has been 
applied to WMAP-SDSS [69] and Planck-SDSS (Li, Zhang & Jing, submitted to ApJ) to 
search for the dark flow induced kSZ effect. In addition to correlations with the distribution 
of galaxies, one also expects a tight correlation between the kSZ effect in the CMB and the 
lensing field (e.g. [87]) or other tracers of large scale structure such as the 21cm intensity 
maps. For brevity, we focus in this work on cross correlating with galaxies. 

For galaxies in a redshift bin [ 2 ;i, 2 ; 2 ], the surface overdensity is dT, = figdgdz. Here, 
fig is normalized such that figdz = 1. We consider a typical bin size of Az = Z 2 — zi ~ 0.2. 
For redshift bins dehned by photometric redshift, we shall take the broadening of redshift 
width by photo-z error into account. Ves{n) only varies slowly with sky direction. So to 
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Figure 11. The kSZ power spectrum for B = 10 Solid, dot, shot dash and long dash lines have 
9f, = 0°,45°, 135°, 180° respectively. Notice that when 9^ = 90°, the kSZ effect vanishes. 


excellent approximation, over degree scales, it can be treated as constant. Hence for a 
degree-sized patch centered on the direction fip, the cross correlation of galaxies with kSZ 
can be well approximated by [67] 


I^Crg(^lhp) =-^ y Alg{k, z)ng{z)xiz)vesinp, z)dz . (3.15) 

Here the electron density-galaxy density power spectrum (variance) A^g{k, z) is evaluated at 
k = following the Limber approximation. We can dehne a temperature fluctuation 

AT = .^^C'TgTcMB/(2'/r). We then have 


AT = -0.167/ik X 



Ves{ne,Ze) 

lOkm/s 


(1 -|- Zg) Xe^gi^Zg'jdz^ . 


(3.16) 


Notice that fig is normalized such that ng{z)dz = 1. The cross correlation depends on 
not only the bubble collision parameters, but also on the galaxy redshift and galaxy bias. 
For brevity, we only show the result for case B with 9^ = 0 and B = 10“^, for five redshift 
bins of [0.01,0.3], [0.9,1.1], [1.9,2.1], [2.9, 3.1] Sz [3.9,4.1] (Fig. 13). Results for other values 
of B scale with H/10“^. We adopt a simple linear bias for galaxy overdensity, 6g = bgSm 
and adopt hg = 1 unless otherwise specified. AT dehned here is not only proportional to the 
kSZ signal in the given redshift range, but also proportional to the density fluctuation in the 
given redshift range. Hence, comparing AT over different redshifts can be misleading. It is 
better to interpret this as the kSZ signal weighted by the r.m.s fluctuation of galaxies in the 
given redshift bin. 
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Figure 12. A schematic of the correlation between kSZ and the distribution of galaxies. In each 
redshift bin (vertical), the kSZ effect maps the density of free electrons i5e(0, ^e) contained in galaxies 

(right panel) to CMB temperature (left panel), with a varying amplitude determined by the collision- 
induced velocity field Vesid, </>, Ze) (e.g. as shown in Figs. 3-6). The measured kSZ effect in CMB alone 
would be a sum over all redshift bins, however the kSZ-galaxy correlation allows one to reconstruct 
the velocity field as a function of redshift; this is kSZ tomography. 


4 Using kSZ tomography to search for bubble collisions 

The kSZ-galaxy cross correlation has unique directional dependence, due to the directional 
dependence in Ueff- This characteristic directional dependence is the key to isolating the weak 
kSZ effect induced by a bubble collision from overwhelming contaminations such as primary 
CMB and cosmic infrared background. This section focuses on understanding this directional 
dependence (§4.1), and utilizing the directional dependence to search for bubble collisions 
(§4.2). 

4.1 The directional dependence of kSZ-galaxy cross correlation 

To detect a characteristic directional dependence of the kSZ-galaxy cross correlation, in 
observations we can split the survey sky into patches of a few degrees in size. The measured 
cross power spectrum depends not only on the multipole but also the relative direction 
between the patch center and the bubble collision center. We denote Ci{t} as the measured 
cross power spectrum in the i-th patch centered on h*. For redshift bins of size Az ~ 0.2, 
Vef[{n,z) does not vary significantly over the patch and over the redshift range. We can 
therefore write Ci{i) as a spatially modulated power 


Ci{£) = C'pivot(^|hpivot)/(hi) . 


(4.1) 






10 100 1000 lo-* 10 100 1000 10< 


Multipole 1 Mullipole 1 

Figure 13. The kSZ-galaxy cross power spectra for B = 10“^ toward the center of bubble collision 
{Of. = 0°). Solid, dot, short dash, long dash and dot-dash curves have galaxies in the range [0.01,0.3], 
[0.9,1.1], [1.9,2.1], [2.9,3.1], [3.9,4.1], respectively. 


where we have chosen to normalize the power against a pivot direction fipivot) nominally 
chosen to be the collision centre. Since the velocity as a function of direction does not vary 
strongly in a degree-size patch, we can neglect the i dependence in / that comes with the 
small angular scale variation of the velocity. From Eq. 3.15, we can then approximate / as 


h 


f{ni) ^ 


dzneaugxvesihi, z) 
dznearigXVefiinpimt, z) 


(4.2) 


In Fig. 14 we show / computed from the definition Eq. 4.1; as can be seen in this figure, 
neglecting the i dependence in / is an excellent approximation. Furthermore, it shows that 
Eq. 4.2 excellently describes the directional dependence. Therefore, given the parameters 
A, B, he and Xc, we can predict the variation of cross power from one patch to another, 
with no need of detailed modelling of electron-galaxy cross correlation. In particular we have 
f{9e = 7r/2) = 0, due to the azimuthal symmetry. Later on we will use Eq. 4.2 for designing 
optimal weighting to boost the cross-correlation signal. 


4.2 The optimal weighting for bubble collisions 

For a fixed redshift bin, we combine all available patches on the sky to obtain a linearly 
weighted cross correlation 

C^{£) = J2cf'°^{e)Wi . (4.3) 

i 

Here the measured cross power spectrum is C'?*’® = Ci + Ni where Ni is the sum of both 
statistical and systematic noise. The optimal weighting depends on the priors we place on 
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Figure 14. The directional dependence of kSZ-galaxy cross correlation f{9e\£). The panels are for 
[0.01,0.3] (bottom left), , [0.9,1.1] (bottom right), [1.9, 2.1] (top left) and [3.9,4.1] (top right). In each 
panel, there are four curves for £ = 1000 (solid), 3000 (dot), 9000 (short dash) and the approximation 
Eq. 4.2 (long dash), respectively. All the curves with A = 0 and B ^ 0 are normalized to unity at 
9e = 0°. For clarity, curves with A ^ 0 and B = 0 are normalized to — 1 at = 0°. One exception is 
the redshift bin [0.9,1.1]. It has AT{9e = 0°) ~ 0. So we normalize with the value at 0e = 45° instead. 
Curves of different £ overlap over each other and are barely distinguishable against each other. This 
shows that (1) f{9e\£) has negligible £ dependence {f{9e\£) — f{9e)) and (2) the approximation Eq. 
4.2 describes f{9e) to excellent accuracy. 


the noise term. We discuss two limiting cases. One is the optimistic case that all large scale 
structure related components (other than kSZ) in the CMB maps have been removed without 
bias by combining multiple frequency bands. In this limit, the expectation value of W is zero 
((W) = 0)- The r.m.s error of the weighted correlation is 

= . (4.4) 

i 

Here, = {Nf). The normalization of W is irrelevant. However it is convenient to fix it in 
order to find the optimal weighting. We choose the normalization 

{£) = Cpi™t(^) ^ ITJi = 1 . (4.5) 


The optimal weighting is then 


Wi = 


hl<y\ 


E filer: 


(4.6) 


The second case is more conservative, and perhaps more realistic. In this case, there are 
residual large scale structure related components other than kSZ in the CMB maps. Therefore 
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(Ni) 7 ^ 0. With a modest prior that these systematic contaminants are statistically isotropic, 
we need VFj = 0 to eliminate these contaminants in the cross correlation. The optimal 
unbiased weighting must minimize 




(4.7) 


Here, af = (Nf) — (W)^- are the two Lagrangian multipliers for the constraints 

Yli ^ifi = 1 and Wi = 0 respectively. The solution is 


[h - {f)s]/af 

* ms - (/)|] E. ■ 


(4.8) 


{f)s is the average over the survey area with the weighting as follows. 


{f)s 


■ 


(4.9) 


The optimal weighting derived earlier for the case of dark flow [67, 69] is a special case 
corresponding to identical noise of different patches. The weighted noise is 




1 1 

ms - (/)!"" e;w ■ 


(4.10) 


The ratio of two weighted noises (case one versus case two) is 1 — {f)'sl {f‘^)s < 1- This 
inequality confirms our expectation that case two is more conservative and hence has larger 
error. 

To be conservative, we will consider the second case in our forecast. The total signal to 
noise combining measurements at all £ in a given redshift bin is 


N 


2 

tot,2: 


E^pivot(^) [{f)s-{f)l] 




(4.11) 


This assumes that errors in different i bins are uncorrelated. Finally we need to combine 
measurements for all redshift bins to obtain the final forecast for signal to noise 


N 


tot 


E 


N 


- tot.z 


We adopt the Gaussian approximation to evaluate Uj, 

CT9,i + ^ cf“)(Cg + Cf) 


a? = 


2M£/sky,i 

(CCMB + cf ^ cm(Cg + Cf) 

2M^/sky,i 


(4.12) 


(4.13) 


Here, B = B{£) is the beam of CMB experiment. is the instrument noise power spectrum 
at the i-th patch of the sky. Due to variation in the exposure, varies across the sky. 
is the power spectrum of the total kSZ effect whose major contribution is the conventional 
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Figure 15. The forecasted constrain on bubble collision. For the future CMB experiment, we adopt 
a FWHM=2 and a noise level of 1/iK per square arcminute. The overlapping sky area is taken as 
/sky = 0.5. Points with filled circles are for the case ot 0 and B = 0. Points with open square are 
for the case of A = 0 and B ^ 0. For the angular resolution of Planck, primary CMB overwhelms. 
So Planck in combination with LSST has limited capability to detect the kSZ generated by bubble 
collision. A future CMB experiment with better resolution can improve the measurement by orders of 
magnitude. We estimate that a PRISM-like experiment has the capability of measure A and B with 
r.m.s error of 10“®. Constraints on B are systematically better than constraints on A when Zc < 4, 
simply reflecting the fact that when A = B, the case B has on the average larger curvature potential 
inside of our horizon. The two constraints at = 4 are roughly the same, for the reason that they 
have roughly the same mean curvature perturbation inside of our horizon. 


kSZ effect including patchy reionization. is the power spectrum of residual foregrounds, 
which includes all large scale structure related contaminants such as CIB, the relativistic ther¬ 
mal SZ effect and the cosmic radio background. Since this relies heavily on multi frequency 
information, whose effectiveness depends on exposure time, this may also vary with direction. 
Variation of the gray-body power index and the temperature of CIB over the sky can add 
more directional variation. Given uncertainties in the modelling of the conventional kSZ effect 
and foregrounds, we will adopt a flat spectrum /^(C''^^^-|-C'^“)/(27r)(2.73K)^ = lO^K^. Cg is 
the galaxy auto power spectrum in the given patch of the sky and the given redshift bin. 
is the associated (statistical) noise. We only consider shot noise so that = 47r/sky,i/Ag,i. 
/sky,i is the fractional sky coverage. Ng^i is the total number of galaxies in this patch of 
the sky in the given redshift bin. The last approximation of Eq. 4.13 holds since the kSZ 
effect induced by bubble collision is always subdominant/negligible to sum of other CMB 
components / contaminations. 

Given that the bubble collision induced kSZ signal is weak, we require both advanced 
CMB experiments and galaxy surveys in order to be sensitive to bubble collisions. Fig. 15 
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Figure 16. The relation between Xc and Zc (left) and 9c and Zc (right) in ACDM over the range 
in redshifts in Fig. 15. Xc is the x-coordinate of the location of the causal boundary and Zc is the 
corresponding redshift [xc = xi^c))- 9c is the angular radius of the bubble collision on the primary 
CMB sky {cos 9c = Xc/xdec)- 


shows the forecasted constraints on A {B), for the combination of Planck and LSST [88]. 
For LSST, we adopt fig = 30/arcmin^. The redshift distribution is adopted as ng{z) oc 
2 ;^ exp(—z/z*) and 2 ;* = 0.5. We assume that LSST and Planck overlap over 50% of the sky. 
This combination has the capability of detecting the bubble collision if max[A,B\ > 10“^. 
Planck has a limited angular resolution of ~ 10 (£ > 2000), so it misses the majority of the 
kSZ signal generated by the bubble collision, which peaks at ^ ~ 10^. To capture the majority 
of it, the angular resolution of CMB experiment should be improved to ~ 2 . Therefore we 
consider a PRISM-like CMB experiment [89] with FWHM=2 and a low noise level of 1/rK^ 
per arcmin^. Combined with the LSST survey, the constraint on A and B can be improved 
by one order of magnitude to (Ja/b ^ 10“^. 

Two major uncertainties in this forecast, besides the bubble collision parameters, are 
the level of residual foreground contamination and the level of non-Gaussianity of residual 
foreground near the angular resolution limit. The Hrst can be alleviated by cleaner foreground 
removal with more frequency bands and lower instrumental noise. A lower residual foreground 
will make the measurement more noise dominated at small scales, and therefore alleviate 
uncertainties caused by non-Gaussianities of residual foreground. 

4.3 Comparison with constraints from the primary CMB 

We can compare this sensitivity to what has been achieved using CMB temperature [48, 
51] and what has been forecast for CMB temperature and polarization [23]. These works 
all analyzed only the linear template, so only provide constraints and projections for the 
amplitude A. CMB studies work with the angular radius of the collision on the CMB sky, 
given by 

cos 9c = . (4.14) 

Xdec 

where Xc is the comoving position of the causal boundary and Xdec — is the comoving 

distance to the surface of last scattering. To assist in comparison between observables in 
Fig. 16 we show the relation between Xc, Zc, and 9c for the Planck best-fit A CDM cosmology. 

The sensitivity of the hierarchical Bayesian method of Ref. [48] is given by the 50% 
completeness curve for the detection of candidate bubble collisions (Figure 5 of Ref. [48]). 
The constraint is expressed in terms of zq, the central amplitude of the collision in the 
CMB temperature. In the Sachs-Wolfe approximation, this is one third of the primordial 
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Newtonian potential i&j, and in the notation of Ref. [48] we can identify A ~ Szq. Over 
the range of 9c relevant for Fig. 15, the completeness curve is fairly flat, and we obtain a 
sensitivity of j^j ~ 1.5 x lO”'^. The analysis of candidate features in WMAP7 data failed to 
provide evidence for bubble collisions, we can interpret this sensitivity as a constraint. 

Ref. [51] performed an analysis of the WMAP7 data using different methods, obtaining 
the constraint —2.10 x 10“^ < aHQ^ {sm 9c)‘^^^ < 2.13 x 10“^ where a is the amplitude of 
the linear template in the comoving gauge. On super horizon scales, we can relate this to A 
in the Newtonian gauge by a factor of 3/5. We then obtain a constraint on the amplitude 
A of |A| < 1.2(sin x 10““^. The angular dependence is not important over the range 

of 9c probed by kSZ tomography, so we can express the constraint as |A| < 1.2 x 10“^, in 
agreement with Ref. [48]. 

Extra information is carried by CMB polarization [23]. Collisions on the angular scales 
probed by kSZ tomography could be detected at similar signal to noise in CMB polariza¬ 
tion as in CMB temperature. The cosmic-variance limited CMB signal to noise using both 
temperature and polarization will therefore be close to the constraints quoted above. 


4.4 Constraints on the lagrangian underlying eternal inflation 

Auxiliary information is needed to interpret our constraints as constraints on fundamental 
parameters in the scalar held lagrangian underlying eternal inflation, as outlined in Sec. 2. 
The most important auxiliary observables are the energy density in curvature and the 
scalar to tensor ratio The fundamental limit on observing curvature is Qk ~ 10“^, 

when it becomes indistinguishable from a long-wavelength curvature perturbation. A next- 
generation CMB polarization experiment can in principle reach sensitivities for the scalar to 
tensor ratio of r ~ 10“^ [90], while futuristic space-based gravitational wave detectors may 
probe down to r ~ 10“® [91]. 

The first case to consider is that of collisions between identical bubbles, in which case 
we can identify B = 214^^® (1 — cos Axsep)^ /15. In this scenario, kSZ tomography can probe 
the entire observable range of 14/^®. If curvature is detected in the future, kSZ tomography 
can completely rule out the presence of identical bubble collisions in our observable universe. 

For collisions between non-identical bubbles, the story is not as clear-cut. Generally 
speaking, we can only constrain combinations of the fundamental parameters 




(1 - cos AXsep) < 


5 / r-obs , 
2]/ 8I4f® 


(4.15) 


rz. - l^ __ A ... ^2 ^ I \ ^R\ 

V^coll ^obs COsAXsep) < ^ Qohs (4-16) 

where Aiimit and Rumit are limits on these amplitudes obtained through kSZ tomography 
(or another method). The first constraint involving depends on the ratio of robs and 

14^'’®, both of which could be unobservably small. Smaller values of robs correspond to small- 
field models of inflation, which arguably need more tuning to get a large number of e—folds, 
and therefore one might not expect an arbitrarily small 14^*^®. This is the only theoretical 
guidance the authors imagine that one might obtain in the absence of a more fundamental 
understanding of inflation. However, in the optimistic scenario where both curvature and 
tensors are detected, meaningful limits on fundamental parameters can be set using kSZ 
tomography. 
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5 Discussion and Summary 


In this paper we have examined how collisions between bubble universes, a possible relic of 
eternal inflation, might be observed using the kinetic Sunaev Zel’dovich effect. The kSZ effect 
we consider is a census of the CMB dipole observed by free electrons on our past light cone, 
which inhabit nearly every vantage point in our observable universe. Because the kSZ effect 
is a census, it has the power to probe large-scale inhomogeneities far better than the cosmic 
variance limited results associated with observations made from one position. Unfortunately, 
being a secondary effect, the amplitude of the kSZ signature expected in the CMB is tiny, of 
order 1/xK. However, a crucial observation is that free electrons trace large scale structure, 
giving rise to a strong correlation between the kSZ effect and large scale structure in the 
universe. This redshift-dependent cross-correlation, known as kSZ tomography, can be used 
as a sensitive probe of ultra-large scale inhomogeneities. 

The effects of a bubble collision are encoded in a planar-symmetric curvature pertur¬ 
bation, and for the class of observers in the vicinity of the causal boundary of the collision, 
can be described by three constants [25, 29]: the position of the causal boundary (comoving 
position Xc or equivalently the redshift Zc), the amplitude of a linear potential perturbation 
(A), and the amplitude of a quadratic potential perturbation {B). The angular dependence 
of the kSZ effect produced by bubble collisions depends on the parameters A, B, and Zc- This 
makes the kSZ effect a powerful discriminator between models, and is essential for isolating 
the kSZ signature from contaminants. 

The main result of our analysis are summarized in Fig. 15. We forecast that Planck and 
LSST can put useful constraints ~ 10“^ on the bubble collision parameters A and B 

over the allowed range of Zc- The primary source of this limit come from the angular resolution 
and instrumental noise of Planck. For LSST and a next-generation CMB experiment such as 
PRISM, the sensitivity can be improved by roughly an order of magnitude to (Ta/b ~ 10“^. 
The in-principle limit on sensitivity, in the absence of instrumental noise and assuming perfect 
foreground subtraction, can be estimated as roughly aA/B ~ ^ Nearly cosmic variance 

limited constraints exist for A from the primary CMB [46-49], yielding A < 10“'^. Therefore, 
kSZ can in principle beat the cosmic-variance limited constraints from the primary CMB by 
nearly 6 orders of magnitude. To put this in perspective, using e.g. A g -^extra ^ 

where Nextra is the number of extra e-folds beyond the number necessary to solve the horizon 
problem, it would be possible at this sensitivity to probe the initial conditions for inflation 
up to 23 e—folds beyond what is necessary to solve the horizon problem. 

To connect the constraints on A and B to constraints on the fundamental parameters in 
the scalar field lagrangian underlying eternal inflation, ancillary data is needed. In particular, 
a detection of curvature and primordial tensor modes would allow for a direct constraint on 

^This number is obtained with several approximations and can only be treated as an order of magnitude 
estimation. (1) We neglect primary CMB. This is indeed valid at ^ > 4000 where other components dominate. 
(2) we assume perfect foreground removal and vanishing instrument noise. (3) we assume sufficiently high 
galaxy number density such that <C Cg. Therefore we can neglect all terms except and Cg in 

Eq. 4.13. (4) we assume a perfect cross correlation between the kSZ effect in a given redshift bin and the 
galaxy overdensity in that redshift bin, namely the cross correlation coefficient r = 1. These approximations 
allow us to sum over all £ and 2 : bins and obtain (S/N)'^^^ ^ ^max/skyC'b®bbie/C''^®^. We first consider case 
B. Using f^C^®gbie/(2vr) ~ B x lO^K (Fig. 9) and j(2 -k) - l/rK for all we obtain (5/N)L 

10^°B(^niax/10'*)^/sky, Valid for |B| <C 1. For £max = 10^, we obtain a sensitivity of ob ~ 10“^°. Constraint 
on A is of the same order of magnitude. The reason is that the kSZ power spectrum for case A is comparable 
to case B when A = B (Fig. 10 versus Fig. 11), and it scales linearly with A. Nevertheless, constraint on A 
can be a factor of 2 weaker than B due to a factor of 2 weaker signal (Fig. 10 versus Fig. 11). 
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the width of the potential barrier describing the collision bubble and the relative energy scale 
of inflation inside the collision and observation bubbles. In the most constraining scenario 
~ 10“^, robs ~ the in-principle sensitivity of the kSZ effect would allow for 

constraints on the properties of the collision bubble down to scales of order 10-10^ TeV 
< 10^ TeV, < 10 TeV). These scales are potentially relevant 

for supersymmetric or grand unified theories which give rise to eternal inflation. 

Bubble collisions in eternal inflation are a predictive theory of inhomogeneous initial 
conditions for our Universe. In this respect, they provide a testbed for addressing more 
generally what we might learn about the observability of inflationary initial conditions by 
providing a benchmark model. Clearly, kSZ tomography has enormous potential to probe 
the initial conditions for inflation. 

Acknowledgment. —PJZ was supported by the National Science Foundation of China 
(Grant No. 11025316, 11320101002, 11433001), National Basic Research Program of China 
(973 Program 2015CB857001), the Strategic Priority Research Program ’’The Emergence of 
Cosmological Structures” of the Chinese Academy of Sciences (Grant No. XDB09000000), the 
NAOC-Templeton beyond the horizon program, and the key laboratory grant from the Office 
of Science and Technology, Shanghai Municipal Government (No. 11DZ2260700). Research 
at Perimeter Institute is supported by the Government of Ganada through Industry Ganada 
and by the Province of Ontario through the Ministry of Research and Innovation. MCJ is 
supported by the National Science and Engineering Research Gouncil through a Discovery 
grant. PJZ and MGJ thank the support of National Science Eoundation Grant No. PHYS- 
1066293, the Simons foundation, and the hospitality of the Aspen Center for Physics where 
this work was initiated. 


A Linear Evolution of gravitational potential and velocity 


The bubble collision generated perturbation is a mixture of both superhorizon and subhorizon 
modes. Therefore the exact evolution of T and v has to be solved numerically. We adopt 
the conformal Newtonian gauge and solve the following equations for radiation and cold dark 
matter in a flat universe [92], 


dQr,o 

dx 

dQr,i 

dx 

d5 

dx 

dv 

dx 

dT 

dx 


k ^ ^ 

aH ~ dx 
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0r,O — 


3aH ’ 3aH 
. k dd' 

1^ 

2 772 




(DmO + iklra ^0r,o) - T - 


k^ 


3a? 




(A.l) 


We have changed the argument from the conformal time r to x = In a. We numerically solve 
the above five equations for the variables ( 0 ^, 0 ; ©r,!) d, —fu, <h) using a modified version of 
the CMBEAST package [93]. The adiabatic initial condition is 


1 3 

0r,o = 2^0 = 0, u = 0 


(A.2) 
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(x-xJ/rj, 

Figure 17. The velocity field at ^ = 0, for case A and case B. The data points are the numerically 
evaluated velocity field and the solid lines are the analytical prediction assuming superhorizon evo¬ 
lution. The two agree excellently. Hence we will adopt the superhorizon evolution throughout the 
paper, which allow us to address many calculations analytically. 


We find that, although there are subhorizon modes in the initial perturbation, the overall 
evolution of $ and v are excellently described by the superhorizon evolution (Fig. 17). For 
superhorizon modes, the evolution in potential at early time where the cosmological constant 
is negligible is [79, 92] 


,,, , , ,,, le^T+y + 9y3 + 2y2 - 8y - 16 

^'SH(a) = - 


(A.3) 


Here, y = a/ueq and Ueq = ^R/^m is the scale factor at the epoch of the radiation-matter 
equality. When y S> 1, 'hgH —^ 'I'sH,i9/10. On the other hand, the evolution of gravitational 
potential at late time where radiation is negligible is 


5 E{a) 

^'(a) oc -Qrn - 

2 a 


f 


da 


E‘^{a)a? 


(A.4) 


Here E{a) = sj is the normalized Hubble parameter and by setting = 0. 
This result holds for both the subhorizon and superhorizon modes [94]. The integral above 
can be excellently fitted by [95] 


9{a) 


_5H^(a)/2_ 

- QAia) + (1 + ^^)(l + 


(A.5) 
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Hence the evolution of a superhorizon mode over all cosmic epoch can be well approximated 
by 


Dx^[a) — 


^'SH(a) 16V1 + y + 9y^ + 2y^ - 8 y - 16 


^SH,i 


lOy^ 

161/1 + y + 9y^ + 2y‘^ - 8 y - 16 
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— 
1 a 


f 


da 


E^a)a^ 


lOy^ 
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On the other hand, we have [79] 
2a^c^H{a) y 


V = —- 


4 + 32/ 


V^ + 


din a 


2a^c^H{a) y 
HlVLm 4 + 3y 


D\^ + 


dD\^ 
din a 


We define the velocity growth rate function as 

2a^H{a) y 


D^{a) = 


4 + 3y 


Zl\|/ + 


dD^ 


din a 


and Dq, are related by 


= 


1 r Il^(a) 

a .in aEJ 


da 


(A. 6 ) 


■ (A.7) 


(A.8) 


(A.9) 


We find that Eq. A .6 and A.7 describe the evolution of the given gravitational potential 
(Eq. 2.2) and the induced velocity excellently (Fig. 17). Numerically we find that subhori¬ 
zon modes slightly smooth the potential and velocity near the causal boundary. But these 
modes are negligible for the purpose of this paper. Therefore we will adopt the above an¬ 
alytical results for superhorizon modes to follow the evolution of bubble collision generated 
perturbation. 
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